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Abstract. We present the results of a 0.86 square degree CCD photometric survey of the open cluster NGC 2516, 

^ ■ which has an age of about 150 Myr and may have a much lower metallicity than the similarly- aged Pleiades. Our 

f~--. ' BVIc survey of cluster members is complete to 1/ ~ 20 and is used to select a preliminary catalogue of 1254 

On . low mass (0.2 < M < 2.0Mq) cluster candidates, of which about 70-80 percent are expected to be genuine. 

^^ ' After applying corrections for contamination by non-members and adding data for higher mass stars from the 

t^^ I literature, we investigate the cluster binarity, luminosity and mass function, mass segregation and total mass. We 

^— ^ ■ find a binary fraction of 26 ± 5 percent, for A to M-type systems with mass ratios between 0.6 and 1, which is 

very similar to the Pleiades. The mass function is metallicity and evolutionary-model dependent, but consistent 

.<^ , with a Salpeter-like law {dN/dlogM oc M"°, a = -1-1.47 ± 0.11 or a = -1-1.67 ± 0.11 for the solar and half-solar 

f| ■ metallicity models of Siess et al. (2000), and a = -1-1.58 ± 0.10 for the solar metallicity models of D'Antona & 

Ph' MazziteUi (1997)), for 0.7 < M < S.OMq. At lower masses (0.3 < M < 0.7Mq) there is a sharp fall in the mass 

function, with a — —0.75 ±0.20 or a = —0.49 ±0.13 (for the solar and half-solar metallicity models of Siess et al.), 

5_j ' and a = —1.00 ± 0.18 (for the solar metallicity models of D'Antona & MazziteUi). The true stellar mass function 

"ti I might have a values up to 0.4 larger if account were taken of low mass stars in unresolved binary systems with 

j^ . mass ratios less than 0.6. The falling mass function of NGC 2516 at lower masses seems inconsistent with the 

much flatter mass functions derived from comparable data in the Pleiades and field populations. This deficit of 

lower mass, fainter stars is also seen in the observed luminosity function. We rule out incompleteness as the cause 

^^ ■ of this discrepancy, but demonstrate that mass segregation is clearly present in NGC 2516, with more than half 

Vh ' the low-mass (< 0.6Mq) stars likely to lie outside our survey area, but the vast majority of high-mass (> 1.5Mq) 

stars included. Taking this into account, it is probable that the whole-cluster mass functions for NGC 2516 and 

the Pleiades are similar down to 0.3Mq. The mass of NGC 2516 stars with M > O.3M0 inside our survey is 

950 — 1200Mq, depending on metallicity and what corrections are applied for unresolved binarity. Correcting for 

mass segregation increases this to ~ 1240 — 1560Mq, about twice the total mass of the Pleiades. If NGC 2516 

and the Pleiades do have similar mass functions, then less massive stars and brown dwarfs contribute about a 

further 15 percent to the mass of NGC 2516 and we predict a cluster population of about 360-440 brown dwarfs 

with 0.030 <M < O.O75M0. 

Key words, open clusters and associations: individual; NGC 2516 - stars: mass function, luminosity function - 
binaries: general 

1. Introduction light element abundances in cool stars with convective en- 

velopes. The ability to study samples with common, well 
Observations of co-eval stars in open clusters play a vi- determined, ages, distances and compositions is the key to 
tal role in investigating the low-mass stellar initial mass ^^^^^ usefulness. A careful census and identification of clus- 
function and defining the physical processes that drive the ^^^ members is usually a prerequisite of all such studies: 
evolution of rotation, magnetic activity and photospheric ^^^ j^ ^^^^^ ^^ p^^^^^^ contamination of sample proper- 

~Q J a ~Z , , D -r, T fl- • ties by interloping non-members that have different prop- 
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Prior to 1990, most work in this area focussed on 
the nearby, well studied clusters such as the Hyades and 



Recently, new X-ray studies by the ROSAT high resolu- 



Pleiades, for which there was ample pre-existing photom- 
etry and proper motion information. However, in the last 
decade it has been realised that a wider range of clusters 
need to be studied in detail. The reasons for this are: 

1. To properly sample a wide age range, from clusters just 
emerging from their parent molecular clouds at 10 Myr, 
through mature clusters like the Hyades at 600 Myr to 
the rarer old clusters such as M67 with ages of a few 



tion imager (Micela et al. 2000), Chandra (Harnden et al 



2001) and XMM-Newton (Bciortino et al. 2001) have de 



Gy 



voted considerable amounts of time to observing NGC 
2516. These more recent studies have used the photomet- 
ric catalogue described in this paper to identify cluster 
X-ray sources. 

So far, systematic membership studies in NGC 2516 
have been limited to fairly bright stars or small areas. 
An unbiased, but severely incomplete list of cluster mem- 
bers down to about V = 15 is presented by Jeffriei 



2. Ass uming that the properties of stars in one cluster 



et al. (19971) , based upon positio n in the V vs B — V 



are 



representative of all such stars at a similar age is 



colour- mag nitude diagram (CMD). [Hawley, Tourtellot, fc 



pre carious and neglects the possible influence of com- 
position differences or birth conditions on behaviour 
in later life. 
3. To investigate dynamical effects such as mass segre- 
gation on cluster evolution and the possible universal- 
ity of properties such as binary fractions and initial 
mass functions requires accurate membership studies 
in clusters with a wide range of ages and stellar densi- 
ties. 

In the last few years the rich, southern Galactic open 
cluster, NGC 2516 (RA= 07''58'", Dec= -60'^45'", / = 
274°, b = —16°), h as become a key obj ect in the study 
of low mass stars. [Meynet et al. (1993| ) give an age of 
141 Myr, in comparison with ages for the better studied 
Pleiades and a Per clusters of 100 Myr and 52 Myr re- 
spectively. Photometric studies of early-type cluster mem- 
bers were presented by pox (1955 ), Eggen (1972[ ) and 
Dachs & Kabus (1989), yielding a mean reddening esti- 
mate E(B-V)=0.12 and a distance of about 400 pc. A ma- 
jor X-ray study was undertaken using the ROSAT satellite 
by [Jeffries, Thurston, fc Pye (1997| ). 159 X-ray sources 
were identified within a 20 arcmin radius of the cluster 
centre, 65 of which could be identified as photometric 
members of the cluster with V< 15, the majority being 
cool, coronally active stars. Using these stars, and compar- 
ing with other clusters, Jeffries et al. (1997| ) showed that 
NGC 2516 has a distance of about 390 pc and a,U — B ex- 
cess that is best explained if its metallicity is substantially 
less than the Pleiades. A model dependent metallicity of 
[Fe/H]— —0.32 ± 0.06 was derived, in agreement with an 
earlier study by Cameron (1985). A sub-solar metallicity 
was also found using an independent B — V vs V — I pho- 
tometric technique by Jeffries, James, & Thurston (1998) 



Reid (1999| ) present photometry and low resolution spec- 



and Pinsonneault et al. (2000), but has yet to be confirmed 
spectroscopically. 

It is this possible low metallicity, together with its rich- 
ness that makes NGC 2516 so interesting. There are plau- 
sible reasons why dynamo generated activity, rotational 
spindown and light element depletion could be profoundly 
affected by a low photospheric metallicity and thus dif- 
fering convection zone properties. The compact size and 
numerous cluster members make NGC 2516 an ideal tar- 
get for fibre spectroscopy and for X-ray satellites such as 
XMM and Chandra, which have limited fields of view. 



troscopy in the cluster d own to much fainter limits, though 
over a very small area. Jeffries et al. (1998 ) present high 
resolution spectroscopy of photometrically selected candi- 
dates between 11 < y < 14.5, establishing definite mem- 
bership through radial velocity measurements as well as 
lithium abundances. This work on brighter objects has 
been continued recently by Terndrup et al. (2001 ). 

In this paper we present a catalogue of CCD photom- 
etry and astrometry down to faint magnitudes {V — 20) 
and over a wide area (0.86 square degrees) in NGC 2516. 
This catalogue will be an invaluable tool for selecting un- 
biased samples of F, G, K and M stars for further study, 
interpreting X-ray observations and studying the dynam- 
ical state and mass function of NGC 2516. A preliminary 



version of this work appeared in Thurston (1999). We 
will restrict ourselves in this paper to a description of the 
data, the construction of a photometric catalogue and a 
membership classification on the basis of this photometry. 
Section 2 outlines the collection and reduction of the pho- 
tometric data, and deals with the astrometric calibration, 
catalogue completeness and comparison with the previous 
literature. Section 3 describes a method for selecting can- 
didate cluster members, which will need to be refined as 
more observational data become available. Section 4 uses 
this catalogue to look at the luminosity and mass function 
of NGC 2516 and search for evidence of mass segregation. 
The catalogue itself can be obtained in electronic format 
from the Centre de Donnees astronomiques de Strasbourg. 

2. Observations 

The data used in this paper were collected at the 0.9-m 
telescope of the Cerro Tololo Interamerican Observatory 
(CTIO) on the nights of 8 9 and 12 January 1995 (lu- 
nar phase 0.52-0.88). A CCD photometric survey was 
completed over an approximate 0.86 square degrees, 
centred upon RA = 07''58'"02.8^ Dec= -60''44'"54^ 
(J2000.0), which comprised of 25 overlapping frames of 
13.5x13.5 square arcminutes, taken with a 2048x2048 pixel 
Tektronix CCD at the Cassegrain focus. The centres of 
the fields were separated by 10.5 arcminutes, leading to 
a square survey area, 55'"36* on a side (see Fig. Q). A 
series of short (30s, 15s, 15s) and long (200s, 100s, 100s) 
exposures were taken through a set of Harris BVI fil- 
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Table 1. Nightly photometric coefRcients, with number 
of standards and root mean square discrepancies 
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Fig. 1. An optical picture from the Space Telescope 
Science Institute digitized sky survey showing the region 
of sky included in our CCD photometry. 



ters respectively. In addition to these target exposures, on 
each night a set of high quality twilight flat fields were ob- 
tained through each filter and many Landolt (1992) stan- 
dard fields were observed over the entire range of airmass 
for which the cluster data were taken. 

All of the frames were bias subtracted and flat-fielded 
with standard methods using the Starlink FIGARO package 
( ^hortridgc ct al. 1999 ). The only complications were that 
the CCD was read out using two amplifiers with slightly 
different bias, gain and readout noise characteristics. This 
was dealt with by reducing the two halves separately and 
then correcting for the gain ratio by requiring that the 
flat flelds were continuous across the amplifier boundary. 
The flat-flelding was tested by flat-fielding median stacked 
night sky exposures and found to successfully remove most 
of the structure in the sky at a level better than a few 
tenths of one percent, apart from several CCD cosmetic 
defects (including a number of bad columns) and a narrow 
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of about 100 Landolt (1992) standards, rejecting stars 
greater than 3cr from the fit at each step until the solution 
converged. 

Equations of the form 



V = v + <Pv{B -V)- KvX + Zv , 
B-V ^ (f>Bv{b ~v)~ KbvX + Zbv , 
V-Ic - (f>vi{v -i)- KviX + Zvi , 



(1) 
(2) 
(3) 



were used to transform the instrumental magnitudes (in 
small letters) onto the Johnson V , B — V, Cousins V — Ic 
system, where </> are the colour coefficients, K the extinc- 
tion coefRcients, X the airmass and Z the zero points 
(the magnitude yielding 1 detected photon per second). 
The best fit values for these coefficients on each night are 
given in Table H^, along with the number of Landolt stan- 
dards measured and the rms discrepancy between their 
measured and published magnitudes. Residuals were ex- 
amined as a function of airmass, colour and time with no 
significant trends seen. On each night we observed several 
red standards with 2 < V — Ic < 2.8, so we expect our 
calibration to be well constrained in this region. 

Each target frame was then analysed. The first step 
was to mask out regions that were saturated due to the 
presence of very bright st ars. The routine daofind from 
the DAOPHOT II package (Stetson 1987; Eaton & Privett 



mask was constructed so that these defective pixels were 
not considered in the subsequent reduction. 

2.1. Photometry 

The Starlink photom package ( Eaton et al. 200C| ) was 
used to determine aperture photometry of the standards 
in a 6 arcsec radius (the stellar FWHM for all our data 
was around 1.3-1.7 arcseconds), which contained about 
95-98 percent of the stellar flux. Nightly transformation 
coefficients, extinction and zero points were determined 
using iterative weighted least squares fits to photometry 



(10 ar( [;second) strip around the outside of the field. A 1996]) was then executed using a 5 sigma threshold and 



with shape rejection parameters set to exclude extended 
objects and cosmic ray spikes. The number of sources 
found per field ranged from 100-200 in the short B ex- 
posures to about 3000 in the long / exposures. 

Aperture photometry was performed using the pho- 
tom routines. Both 3 and 6 arcsec radius apertures were 
used. For the brighter stars, where the statistical errors 
were smaller than 0.01 mag (about V — 17), the larger 
aperture results were used and these also calibrate the 
aperture correction from the small to the larger aperture 
for the fainter stars. The sky estimation was taken from 
the mode in a surrounding annulus. At this stage stars 
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Fig. 2. Top: The V ms B — V CMD for our survey. The s olid line shows the ^iess et al. (2000| ) solar metallicity isoch rone 
we use to select members of the cluster (see Sect. 3^). Bottom: Objects selected as cluster members (see Sect. 3.3 ) 
are shown as triangles (where both B ~ V and V ~ Ic were available) or squares (where only B ^ V was available - 
for just 5 objects). 
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Fig. 3. Top: The V vsV — Ic CMD for our survey. The s olid Hue shows the ^iess et al. (2000| ) solar metallicity isochrone 
we use to select members of the cluster (see Sect. 3^). Bottom: Objects selected as cluster members (see Sect. 3.3 ) 
are shown as triangles (where both B ~ V and V — Ic were available) or squares (where only V ~ Ic was available - 
for 368 objects). 
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Fig. 4. Top: The V — I^ vs B ~ V colour-colour diagram for our survey. Bottom: The 881 objects selected as cl uster 
members and which have both B — V and V — L measurements are shown as triangles. The solid line shows the BiesE 



et al. (2000 ) solar metallicity isochrone we use to select these members (see Sect. p.3[ ). 
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were rejected if masked pixels fell within the object aper- 
ture and if they were close to the frame edges. 

The next stage combined the separate filter measure- 
ments in each field, treating the short and long expo- 
sure sequences separately. The instrumental b,v,i mea- 
surements for a particular field were combined by match- 
ing the pixel coordinates and allowing for small transla- 
tions between the fields. The v object list was used as the 
reference frame (i.e. an object must be detected in v and 
at least one other filter to be incorporated in the final 
catalogue). Mean airmasses were calculated for the B,V 
and V, / pairs and the catalogue of measurements for that 
field transformed onto the standard magnitude system. If 
no B — V measurement was available (faint, red objects), 
then we used a mean relationship between B~V and V—Ic 
in order to estimate B — V for the purpose s of transforma- 
tion to a standard V magnitude (see Sect. 3.3). This latter 
procedure should add no more than ±0.01 magnitudes to 
the V error because the colour term, (py, is small. 

A preliminary astrometric solution was found for the 
centroided CCD positions in the long V exposure of each 
field. This was obtained by identifying many stars from the 
Guide Star Catalogue Version 1.1 (Lasker et al. 1990). A 6- 
coefficient fit to these reference stars, resulted in solutions 
good at the level of 1 arcsecond. The x,y pixel coordi- 
nates of objects identified in the long exposures were then 
transformed to RA and Dec using this solution. Of course 
some bright stars {V < 14) were saturated in the long ex- 
posures and their positions must be taken from the short 
V exposure. This was achieved by fitting a linear trans- 
lation between the long and short exposures, using stars 
detected in both frames. 

Objects were then matched between the catalogues 
from the short and long exposures. Where objects appear 
in both, a weighted mean of the photometry was taken. 
There were generally about 200 well measured stars com- 
mon to both the short and long catalogues in each field. 
The magnitudes of these stars were in very close agree- 
ment. The biggest differences we ever found were about 
0.02 mag, which confirmed our preliminary assessment 
at the telescope that these were very good photometric 
nights. 

We followed a similar procedure to deal with stars in 
the overlapping regions between fields. Again, where two 
(or even more) measurements existed, a weighted mean of 
the photometry was taken. Analysis of the discrepancies 
between magnitudes of stars in these overlap regions is 
our primary estimate of the internal accuracy of our pho- 
tometry. We found no evidence for systematic variation in 
the photometric calibrations between nights or between 
fields at a level greater than 0.02 mag. Table || shows the 
internal error estimates (the rms values for the overlap dis- 
crepancies) as a function of magnitude. Beyond 1^ = 20 
the statistical errors in the photometry rise rapidly. 

Possible causes of error are (of course) the statistical 
errors, but also variations in the point spread function, 
and hence the aperture correction, over the CCD field of 
view (especially in frames not precisely in focus). This 



Table 2. An estimate of photometric errors from stars 
measured more than once in different frames 
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0.032 
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0.040 


253 


0.026 
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0.050 


395 


0.041 


318 
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360 
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0.075 
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0.046 


257 


0.050 


401 


3-20 


0.099 


534 


0.052 


65 


0.055 
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would contribute a term (in quadrature) which was the 
same at all magnitudes - which is approximately what we 
see. Although the short and long exposures in each field 
were taken in one observing sequence, some of the over- 
lapping fields are separated by hours or even nights. It is 
therefore quite plausible that some genuine variability also 
contributes to these errors. That the errors in the V mag- 
nitudes are larger than the colour indices, suggests that 
errors in correcting measured magnitudes to the standard 
star aperture values near the edges of the CCD fields are 
the more likely culprit. If this is the case, the errors in 
Table y are likely to be overestimates for the majority 
of stars with V < 18. For V > 18, the statistical errors 
dominate. 



2.2. Astrometry 

The preliminary positions determined from 6-coefficient 
fits to objects in the Guide Star Catalogue (rms typi- 
cally ~ 1 arcsec) were improved upon usi ng pre-release 



data f rom the SuperCOSMOS Sky Survey (|Hambly et al , 
^001a| ). The UK Schmidt Bj plate from field 124 was used 



(plate number J2978, epoch 1977.223). These data con- 
sisted of a catalogue of ^ 1.4 million objects to Bj ~ 22. 
At the time these data were used, the global astrometric 
plate reductions were based on standards from the Tycho- 
AC catalogue ( Urban, Corbin, fc Wycoff 199S ) and had 
typical residuals per standard of ^ 0.3 arcsec in either co- 
ordinate (for more details of the astrometric reductions 



for Su perCOSMOS Sky Survey data, see |Hambly et al 
[2001b ). For each CCD frame, objects in the photometric 



catalogue were matched with objects found on the pho- 
tographic plate. A matching radius of 5 arcsec was used. 
A 6-coefficient linear transformation was then applied to 
the CCD coordinates with an additional cubic radial dis- 
tortion coefficient as a further free parameter. The optical 
axis was assumed to be the centre of the CCD. Parameters 
were adjusted to get the smallest rms when compared with 
tangent-plane positions on the photographic plate (the 
radial distortion term was almost negligible). The typi- 
cal zero-point shifts applied to the preliminary positions 
were about 0.7 arcsec, with final rms values of around 0.3 
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Table 3. The photometric catalogue for NGC 2516, in- 
cluding J2000.0 positions, V, B — V and V—Ic photometry 
where available, flags indicating whether each star passes 
the three membership tests discussed in Sect. 3.3, a flag 
indicating whether the star is a candidate cluster member 
and finally a flag indicating whether a cluster candidate 
is a binary system. The Table is available in electronic 
format by ftp from CDS at cdsarc.u-strasbg.fr. 
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Fig. 5. The apparent V luminosity function for all the 
stars detected in our survey (solid squares) compared with 
the luminosity function for those stars with V — I^ > 2.0. 



arcsecs, which we expect are largely dominated by uncer- 
tainties in the photographic positions. 

The quality of the astrometry has recently been tested 
with a fi-bre-spectroscopy run on the low-mass cluster can- 
didates (see Sect. |3|), using hydra on the CTIO 4-m 
telescope. Excellent results were achieved over a 40 ar- 
cminute diameter field, using brighter cluster members as 
the fiducial acquisition stars and 2 arcsecond diameter fi- 
bres (Jeffries et al. in preparation). 



2.3. The catalogue completeness 

The photometric/astrometric catalogue is given in Table ra 
(available from the Centre de Donnees astronomiques de 
Strasbourg) and consists of 15495 stars with V magni- 
tudes, B — V and V — Ic colours when available, along with 
their J2000.0 positions and flags indicating their member- 
ship and binarity status (see Sect. 3^). The complete V 
vs B — V and V vs V — Ic colour- magnitude diagrams (the 
BV and VI CMDs) are shown in Figs. and y and the 
colour-colour diagram is shown in Fig. y. We refer here- 
after to this as the CTIO catalogue. 



In Fig. H, the apparent V magnitude luminosity func- 
tion of the whole catalogue is shown in 0.25 magnitude 
bins. There is a clear turnover in the apparent luminosity 
function beyond V = 19.25. Catalogue completeness has 
been investigated by adding simulated stars to our ob- 
served frames and then searching for these stars using the 
same daofind parameters. We find that there are three 
possible causes for the apparent turnover in the luminosity 
function. 

1. Clearly our search threshold imposes a magnitude limit 
on the catalogue. Tests on the long V frames near the 
centre of the cluster reveal that the level of incomplete- 
ness caused by this and also by faint stars being close 
to bright stars is approximately 1.5 percent at T^ = 19 
and 10.7 percent at V^ = 20. This is not sufficient to 
account for the observed turnover. 

2. We require that in addition to being detected in the V 
observation, a star is also detected in at least one of the 
B or / observations. The individual 5-sigma detection 
limits for our long exposures average about V = 20.5, 
/c = 19.3 and B — 20.4 (depending on seeing and 
moon phase). Thus in practice their are only a few 
(faint blue) stars detected in V and B only and the 
vast majority of the catalogue has a, V — Ic measure- 
ment. However, there will be a significant number of 
faint, intermediate colour stars (0.5 < V — Ic < 1.5) 
which are detected in V but not /. Our simulations 
show that for V — Ic = 0.9 (where the bulk of the pop- 
ulation might be expected to lie), 10 and 50 percent 
incompleteness occur at about / = 18.9 and / = 19.2 
and therefore V = 19.8 and y = 20.1. 

3. This latter effect is confused by the probability that 
the turnover in the luminosity function is real. At y — 
Ic = 0.9, the major population of the CMD at y = 20 
will be K-dwarfs, because the cluster has a galactic 
latitude of -15.8 degrees. At 1/ = 20 a typical K-dwarf 
would be at a distance of 7.5 kpc and therefore more 
than 2 kpc below the galactic plane. The scale height 
of such stars is actually less than 1 kpc, so we would 
naturally expect to see a decrease in the luminosity 
function beyond V ~ 18.8. 

The latter two effects become much weaker for redder 
stars. These are invariably detected at / and because M- 
dwarfs are much less luminous than K-dwarfs, the scale 
height argument is also less applicable. The completeness 
limit for faint, red stars in NGC 2516 will therefore be 
determined solely by the first factor. This is illustrated 
in Fig. ra, where we also show the apparent V luminosity 
function for those stars with V — Ic > 2.0. The final out- 
come of this discussion is that we believe our catalogue 
is almost (> 98 percent) complete to F = 19 over the 
whole V vs V ~ Ic CMD. However, we estimate that our 
catalogue is still 90 percent complete to V^ = 20 for faint 
members of NGC 2516 that are selected on the basis of 
V and V — Ic- If we require selection in the V vs B — V 
CMD as well, the completeness level is governed by the 
B frame completeness for red stars and we have 90 per- 
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cent completeness for faint, red NGC 2516 members at 
about V = 18.2, and about F = 19 for the general field 
population at B — V = 0.8. 



2.4. Comparison with otiner photometry 

A number of other photometric studies have been pub- 
lished on NGC 2516, though none that go both as deep 
and cover as large an area as this one. We have compared 
our work with photoelectric {BV) photometry published 
by pilggen (1972 ), photographic and photoelectric {BV) 
photometry published by Dachs fc Kabus (1989) and CCD 
{Vic) photometry published by Hawley et al. (199{| ). 



The |Eggen (1972D and [Dachs fc Kabus (1989D BV pho - 
tometry is discussed extensively by Jeffries et al. (1997D . 
In that paper an optical catalogue is compiled from these 
sources, in which small corrections (of order 0.02 mag.) are 
made to the Eggen photoelectric and the Dachs & Kabus 
photographic values, to put them on a consistent system 
defined by the Dachs fc Kabus photoelectric values. Dachs 
& Kabus estimate 0.02 mag. errors in their photoelectric 
values and Jeffries et al. deduced errors of about 0.06 mag. 
for their photographic values. The final catalogue contains 
568 objects and was estimated to be complete to about 
V = 13.5. We find that there are 500 matches with our 
new CTIO catalogue. The remaining objects either lie out- 
side the area covered by the CTIO survey or were bright 
objects and hence saturated. Hawley et al. (199£ ) present 
a list of 155 candidate members of NGC 2516 that are 
within our survey area. We were able to find 130 matches 
to these within 1.2 arcseconds. The remaining unmatched 
stars are very faint and were unlikely to be detected in 
our survey. 

Graphical comparisons between the catalogues are 
shown in Figs, g, and || where we have divid ed the 



comparison between: the photoelectric data from Eggen 
(1972| ) and [Dachs fc Kabus (1989D (which are given difller- 



ent symbols); the photographic photometry from Dachs 



fc Kabus; and the CCD photometry from Bawley et al 



(1999). The lower panels in these Figs, show the discrep- 



ancies between the CTIO and comparison photometry. In 
all plots, the solid line simply represents equality between 
the CTIO and comparison data, rather than a fit. 

Figure o shows that there is an increasing discrepancy 



towards fainter V magnitudes among the the Dachs fc 



Kabus (1989) and Eggen (1972) photometry in compari- 
son with the CTIO data. This appears to be the case to a 
much severer degree in the photographic Dachs fc Kabus 
data (Fig. M). We are confident that our photometry is 
excellently calibrated in this V magnitude range, with no 
trends at all to T^ = 16. We note that the two labelled 
outliers, DK427 and DK865, are quite close (< 6 arcsecs) 
to companion stars which might easily have made them 
appear brighter in the poorer quality photographic pho- 
tometry. In contrast the B — V data show very good agree- 
ment with our data, although there is some indication that 
the photographic colours are too blue by about 0.05. We 



have inspected our images and can see no problems with 
our photometry of the star called E3 by Eggen (labelled 
in Fig. n). This may well have been misidentified in Eggen 



(1972), which is reinforced by the agreement of our value 



with that obtained by Hawley et al. (1999) for the same 
star. The scatter in the residuals for Figs, o and M are 
in line with our estimates of the errors in these datasets, 
although the Eggen photoelectric photometry appears to 
show more scatter than that of Dachs fc Kabus. 

Turning to Fig. g[we again have reasonable agreement. 
There is a systematic offset of about 0.04 magnitudes in V 
(in the sense that the Hawley et al. photometry is brighter) 
and a suggestion that this reverses for the faintest stars 
in the sample. We note that as Hawley et al. only pub- 
lished data for cluster candidates, all these very faint stars 
are also very red. It therefore suggests that there may be 
minor problems with either their or our colour-dependent 
terms in the transformation equations. We re-iterate that 
we did observe Landolt (1992) stars as red as V — Ic — 2.8 



and we are certainly confident in our calibration to this 
point (and the colour term is very small in any case for 
our CCD and filter combination). This pattern is repeated 
for the V — Ic comparison, where Hawley et al.'s colours 
are blue by about 0.04 magnitudes compared with ours, 
with a definite reversal in the brightest stars. The scatter 
in the residuals of both comparisons is almost precisely 
in accord with the error estimates given by Hawley et al. 
and in this paper. We have no explanation for the very 
discrepant point HTR103 (labelled on the diagram), other 
than perhaps a major stellar flare occurring during Hawley 
et al.'s observation. 



3. Selection of Candidate Cluster Members 

3.1. Selection philosophy 

We have used the photometric catalogue to attempt a pre- 
liminary selection of cluster members based only on photo- 
metric criteria. This selection procedure is especially use- 
ful because it does not rely on any characteristic of the 
cluster members which one might choose to investigate - 
for instance coronal X-ray emission. Such a catalogue of 
cluster candidates will be unbiased with respect to mag- 
netic activity, rotation rate or lithium depletion and there- 
fore provides an ideal starting point for such investigations 



(see 

pool 



Jeffries et al. 1998; Miccla ct al. 200C; Harndcn et al. 



Bciortino et al. 2001) 



We must provide a caveat here; our selection procedure 
is arbitrary to some degree, and will exclude some genuine 
members and will include some non- members. We will 
fashion our photometric selection criteria so as to avoid 
excluding the vast majority of cluster members. The inter- 
ested reader should easily be able to generate membership 
catalogues using their own (possibly more restrictive) cri- 
teria. Ideally these catalogues should then be refined using 
other unbiased indicators of membership such as proper- 
motions or radial velocities. 
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Fig. 6. A comparison of photometry between this paper and photoelectric photometry from Dachs & Kabus (1989) 
[sohd symbols] and Eggen (1972) [open symbols]. The straight line simply indicates equality and is not a fit to the 
data. 



Our CCD data combined with the SuperCOSMOS 
scan of the 1977 epoch Schmidt plate (see Sect. |2^ ) did 
allow a preliminary attempt at proper-motion selection 
for the fraction of the catalogue which had good pho- 
tographic positions. We found that the best accuracies 
achievable were about 6 milli-arcsecondsyear"^ in each 
coordinate. We calibrated the cluster mean proper mo- 
tion using known members from Jeffries et al. (1998 ), and 
found it to be essentially zero within the errors. It was 
soon discovered that even for the most accurate data, we 
were unable to exclude more than about 10 percent of the 
general field background contamination whilst including 
more than 90 percent of the cluster members. We will not 
present these preliminary proper motion results in this 
paper, but will await a second epoch CCD survey which 
should be capable of producing more precise and useful 
results. 

Our selection philosophy is therefore restricted to using 
the two CMDs {V vs B - V and V vs V - /c) to select 



stars close to a cluster isochrone and then to check the 
colour-colour diagram for consistency. 

3.2. Isochrone generation 

We used 150 Myr approximately solar metallicity (z — 
0.02) isochrones generated from the models of piesa 



[Pufour, fc Forestini (200(]| ) over the mass range O.IMq < 
M < 7Mq. These isochrones were provided in the form 
of bolometric luminosity versus effective temperature and 
had to be transformed into the observational plane. We did 
not use the photometric colours provided with the models. 
The critical step here is the relationship between colour 
and effective temperature, which is highly uncertain, es- 
pecially among cooler stars. 

The procedure we adopted is as follows: We obtained 
a catalogue of BVI photometry of the Pleiades (cour- 
tesy of the Open Cluster Database operated by J. Stauffer 
and C. Prosser) and supplemented this with photometry 
[of cooler Pleiades members presented by Stauffer et al 
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Fig. 7. A comparison of photometry between this paper and photographic photometry from Dachs & Kabus (1989) 
The straight hne simply indicates equahty and is not a fit to the data. 



(1998a). V — I indices on the Kron system were converted 
to the Cous ins values using transformations in Bessell & 
Weis (L987). We posit that there is a single colour- Toff 
relationship which applies to all stars of approximately 
the same age as the Pleiades, for which we assume an 
age of 120 Myr, a distance modulus of 5.6, a reddening of 
E{B -V)^ 04, E{V - Jc) = 0. 5 and extinction of 
Av = 0.13 (see [Stauffer et al. 1998a| ; [Stauffer, Schultz, fc 



Kirkpa trick 1998b| ). Note that the exact values of these 



quantities have very little influence on the membership 
selection for NGC 2516 but could affect the deduced dis- 
tance modulus or metallicity for the cluster from isochrone 
fitting. 

An empirical isochrone is fitted (by eye) to the Pleiades 
data (in both BV and VI CMDs), ensuring that the 
isochrone is not biased upward by the presence of binary 
cluster members (see Fig. ||). At about 25 fiducial points 
along these empirical isochrones we determine absolute 
magnitudes and colours, which are then converted to bolo- 
metric luminosities using bolometric correction-colour re- 



lationships from polynomial fits to the empirical data in 
Flower (1996) for B — V, and a c ombination of empiri- 
cal data from Leggctt et al. (1996) for V — Ic > 0.7 a nd 



the atmospheric models of [Bessell, Castelli, fc Plcz (1998| ) 
for V ~ Ic < 0.7. Judging by the scatter around these 
relationships and by comparisons with other bolome tric 
correction-colour relationships (e.g. Monet et al. 1992J ) we 
estimate systematic errors of no more than a few hun- 
dredths of a magnitude (at least over the colour range we 
are interested in). This contrasts markedly with the un- 
certainties if one tries to use a bolometric correction- Toff 
relation, which is often what has been done in the litera- 
ture to convert models to observables. 



Having obtained luminosity as a function of intrinsic 
colour, we then use the assumed age of the Pleiades to ob- 
tain the Toff value appropriate for a particular luminosity, 
by interpolating along the 120 Myr model isochrone. This 
in turn defines a set of colour-Teff points which can then 
be used to transform any other model isochrone. 
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Fig. 8. A comparison of photometry between this paper and CCD photometry from Hawley et al. (1999). The straight 
hnc simply indicates equality and is not a fit to the data. 



The net result is that we have well calibrated 
isochrones over the colour range defined by the fidu- 
cial Pleiades data (or the mass range of the models - 
whichever is more restrictive). In practice this means 
-0.15 < (B - V)o < 1.7 and 0.4 < {V - Ic)o < 4.0. 
The hotter limit to the V — Ic range is defined by a lack of 

V — Ic photometry for hotter Pleiades stars. For the sake 
of defining an isochrone for membership selection, we have 
extended this data back to (V^ — /c)o = —0.1 by transform- 
ing the hot star B — V dat a into V — Ic using relationships 
defined by Johnson (1966) (bet ween B — V and Johnson 

V — I) and then Bessell (1979| ) (between Johnson V ~ I 
and V - Ic). 

The 150 Myr isochrones were transformed into the BV 
and VI CMDs, and reddening and extinction appl ied for 
NGC 2516 (Ay ^ 0.38, E(B - V) = 0.12 - see iDachs 



V - Ic < 0.8 and lA < V - Ic < 2.2 in the VI CMD, 
where contamination by field stars appears to be small (see 
Sect. 3.4). The results are shown in Figs. || and |^, where a 
distance modulus of 8.10 has been applied to both the BV 
and VI CMDs, with estimated uncertainties of ±0.05. 



This procedure is quite robust to the assumed ages of 
the Pleiades and NGC 2516, because for all but the very 
hottest and coolest stars, changing the age by as much as 
50 Myr makes little difference to the isochrones. It is also 
robust to the choice of evolution model, because we require 
the models to fit the Pleiades at a similar age. We obtain 
essentially an identical fi t to the NGC 2516 data using 
150 Myr isochrones from D'Antona fc Mazzitelli (1997 ) 



and with the same distance moduli as above (see Fig. |1C 

Of more consequence is the possibility that t he metal- 
licity of NGC 2516 is sub - solar by a factor of two ([Cameron 



& Kab |us (1989[) and [Jeffries et al. (1997|, |l998[) - and an 1985| ; [Jeffries et al. 1998[ ; [Pinsonncauh et al. 2000| ). If this 



assumed E(y — I) — 0.15). We then adjusted the intrinsic 
distance modulus of the isochrones to match the data, 
particularly for B - V < 0.9 in the BV CMD and for 



is the case then the distance moduli we have obtained 
will be overestimated because lower metallicity stars are 
fainter at the same colour. The effect should be larger 
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Fig. 9. The Pleiades BV and VI CMDs showing how we 
defined the isochrone fits using 3iess et al. (2000| ) models 
and a distance modulus of 5.6 (see text). 





Fig. 10. A comparison of isochrones generated from the 
three models discussed in Sect. 3.2. The solid lines are 



from the solar metallicity model of Siess et al. (2000), the 
dashed lines are the half-solar metallicity models of Biess 
|et al. (200C ) and the dotted lines (virtually indistinguish- 
able from the solid line) arc from the solar metallicity 
models of D'Antona & MazzitcUi (1997). The isochrones 



have been shifted by different distance moduli to match 
the observational data (not shown for clarity - see Figs. || 
and y) . 



in the BV CMD than the VI CMD because the opac- 
ity caused by metal-lines causes blanketing in the B band 
(Alonso et al. 199(;). Pinsonneault et al. (1998) have cal- 
ibrated this effect for F and G stars and indeed used 
it to calculate the metallicities (and mctallicity-corrected 



distan ce mod uli) f or a number of open cluste rs. Jeffries 
et al. |(1998|) and [Pinsonneault et al. (2"00C| ) used the 
same approach to calculate a metallicity for NGC 2516 
of [M/H]=-0.18±0.08 and -0.26 respectively on the basis 
of preliminary BVI photometry of known members (from 
this dataset). 



We have checked the effects of a low metallicity by us- 
ing a 150 Myr isochrone from 3iess et al. (2000| ) with z = 
0.01 (approximately half-solar metallicity, [M/H]'^ —0.3). 
This is again calibrated using the Pleiades photometry. A 
question arises as to whether the colour- Teg relation de- 
rived empirically from the Pleiades is suitable for a lower 
metallicity cluster. The B — V index is sensitive to metal- 
licity for warm stars with partially ionized metal lines and 
is partly the reason that the BV CMD may be changed 
more by metallicity than the VI CMD. However for hot 
stars {B -V < 0.2) and cool stars (i? - V^ > 1.3) this is 



not likely to be the cases (see pastelli 1999| ) . p^eggett et al. 
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(1996) also show that the V — Ic colour index is a good 
temperature indicator with relatively little metallicity sen- 
sitivity for cool stars {V — Ic > 1.5) and line blanketing in 
the V band is not expected to be very important for hot- 
ter stars. Thus to first order, this approach to generating 
a low metallicity observational isochrone should be valid. 
We find that the solar and half-solar metallicity 
isochrones yield comparable fits to the data. As expected, 
the distance moduli required to fit the data are smaller 
for the lower metallicity models. We find distance moduli 
of 7.85 ± 0.05 and 7.90 ± 0.05 for the BV and VI CMDs 
respectively. These distance moduli are in excellent agree- 
ment with the 7.94 ± 0.04 in [Jeffries et al. (1997D and 
7.96 ± 0.17 found by [Jeffries et al. (1998|), but a little 



larger than the 7.77 ±0.10 deduced by Pinsonneault et al. 
(2000|)]and the Hipparcos distance of 7.70 ± 0.16 found 



by Robichon et al. (1999). We emphasize that our errors 
are underestimated because they do not take into account 
uncertainties in the metallicity or reddening. 

The isochrones are comp ared in Fig. pj)[, where we also 
include the solar-metallicity D'Antona & Mazzitelli (1997) 
model. The isochrones have been shifted to the distance 
moduli required to give a reasonable fit to the observa- 
tional data (which is not shown for clarity of compari- 
son between the models). The most important point to 
make is that the shapes of the models are extremely sim- 
ilar. The main discrepancy occurs in the cool part of the 
VI CMD, where low metallicity stars lie just less than 
0.1 magnitudes below solar metallicity stars of the same 
colour. This is fortunate, because it means irrespective of 
which model/metallicity we choose, the selection of clus- 
ter members by photometric means is almost unaffected 
- although cluster properties such as the mass function 
are (see Sect. 12). A more detailed investigation of the 
metallicity, distance and reddening is left to another pa- 
per that uses a more complete sample of spectroscopically 
confirmed NGC 2516 F and G stars ([Terndrup et al. 200l|). 



3.3. Membership criteria 

To select the members we apply the following crite- 
ria, based on the ^iess et al. (2000 ) solar metallicity 
isochrones: 

1. li a. V — Ic colour is available, the star must have V 
within a region bounded by 0.1 mag below the clus- 
ter V vs V — Ic isochrone, and 0.85 mag (to include 
binary stars and allow for the uncertain metallicity - 
see Sect. [3.2[ ) above it. Additionally we allow an extra 
2cr error in V (see Table 0), which also incorporates a 
(usually dominant) contribution from the V — Ic error, 
which is added in quadrature to the V error by assum- 
ing a gradient in the CMD appropriate for a cluster 
member at that colour. 

2. li a B — V colour is available then the star must satisfy 
a similar condition in the V vs B — V CMD. At faint 
magnitudes {V > 17), the cluster isochrone steepens. 
For these stars we additionally allow a cluster member 



to lie within 2aB-v of a band running ±0.05 mag in 
B — V either side of the cluster isochrone. 

3. li B — V and V — Ic are both available, the star must 
lie within ±0.05 mag. (in V — Ic) of the cluster locus in 
the colour-colour diagram (see Fig. Eh. Again, we allow 
an extra slop of 2a which includes contributions from 
both the B — V and V — Ic errors added in quadra- 
ture. The V — Ic vs B — V cluster locus is calculated 
by combining and interpolating the two CMD cluster 
isochrones. Redder than V — Ic > 1.9 this locus steep- 
ens and similarly to test 2 above, we allow a cluster 
member to lie within 2aB-v of a band running ±0.05 
mag in B — V either side of the cluster locus. 

4. Finally, if the star passes test 1, it is classed as a binary 
candidate if it lies more than 0.3 mag in V above the 
V vs V — Ic isochrone. This should select unresolved 
binaries with mass ratios of approximately 0.6-1 (see 
Sect. U). 

The full catalogue contains 15 495 stars. Of these we 
find that 11 114 have a B — V value, 15 310 have aV — Ic 
value and 10 929 have both. 1499 stars pass test 1, 1295 
stars pass test 2 and 920 stars pass them both, of which 
881 pass test 3. However, in order to be classed as a cluster 
member we only require that the star should not fail any 
of these tests. Therefore we include 368 stars which have 
no B — V colour but pass test 1. These are almost all very 
faint, red cluster candidates. We also include 5 stars which 
possess no V — Ic but pass test 2. The final membership 
catalogue therefore contains 1254 candidates. Figures |[, [^ 
and |4| show these cluster candidates with different symbols 
indicating how they have been classified. Of the 1254 clus- 
ter candidates, 403 are possible unresolved binary systems. 
Flags corresponding to the status of stars with respect to 
these tests are appended to the catalogue in Table y. 

3.4. Contamination by non-members 

Whilst we are reasonably confident that our membership 
criteria should have included the vast majority of true 
cluster members (within the bounds of our photometry 
completeness), it is clear looking at Figs. || and [^ that we 
have also included contaminating foreground and back- 
ground sources. As we have taken no offset fields of the 
cluster, we attempt to estimate membership probabilities 
and quantify this contamination by interpolating the dis- 
tribution of cluster non-members in the VI CMD. 

We proceed by removing the cluster candidates from 
the catalogue and then select stars in several strips in the 
VI CMD, both above and below (but parallel to) the fidu- 
cial cluster isochrone (from 1.2 magnitudes below to 2.0 
magnitudes above), but beyond the limits defining mem- 
bership. We fit a smooth function to the density (num- 
ber of stars per V magnitude interval) of contaminating 
stars along vertical strips (colour ranges) in the CMD. 
We achieved the best looking fits using the sum of a con- 
stant density and an exponential decay with an e- folding 
length of 0.5-0.9 mag. An example is shown in Fig. ^ 
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1.68<V-I <2.04 



9 single contaminants 




The results in Table H can now be used to correct sta- 
tistical ensembles for contamination (e.g. in investigating 
the luminosity function - see Sect. 4.1). We caution how- 



Fig. 11. An example of our fitting procedure to determine 
the level of contamination among our candidate members 
(see text). The y-axis shows the number of star per unit 
magnitude found in strips a distance AV (on the x-axis) 
from the cluster isochrone. The example shown is for strips 
bounded by 1.68 < V - I^ < 2.04 (see Tabic |). 



The fitted function is then integrated over the range in 
which cluster members were selected, in order to estimate 
how many contaminating stars we expect to fall within 
the membership selection region. Because the contamina- 
tion tends to decrease quite sharply as we move above 
the cluster isochrone we subdivide this range further in 
an attempt to estimate how much contamination there 
would be among single and binary star candidates sepa- 
rately. This was achieved by integrating the the interpo- 
lating function above and belo w a line 0.3 mag. above the 
cluster isochrone (see Sect. 3^). These figures are then re- 
duced by the number of objects which passed test 1, but 
failed tests 2 or 3, because we argue that these must be 
part of the contaminating sample. Finally we subtract the 
remaining contaminants from the number of cluster can- 
didates in the same magnitude range and divide by the 
number of cluster candidates to estimate the probability 
that a cluster candidate is a genuine cluster member. At 
the same time we can see what advantage has accrued from 
using more than one colour in assessing cluster member- 
ship, by finding what fraction of the contamination in the 
VI CMD has been rejected using B — V. 

The results are given in Table 0, where we have split 
the sample into colour ranges which roughly correspond to 
1 magnitude V intervals for single stars in the cluster. Such 
crude binning is necessary in order to have enough stars 
(especially in the CMD strips above the cluster isochrone) 
to yield reliable fits for the interpolating function. We have 
also assumed that all hotter candidates (roughly brighter 
than V = 11.5) are members. 



see Sect 

ever that these membership probabilities are averages over 
our field of view. Because the cluster is centrally concen- 
trated (see Fig. n^ and Sect, n^ and we expect that the 
background contamination has a uniform spatial distribu- 
tion, membership probability will be higher for a cluster 
candidate close to the cluster centre and lower for a can- 
didate near the edge of our survey area. 

The contamination is worst for 0.87 < V — Ic < 1.33 
(12.5 < V < 15.5) and a glance at Fig. y confirms 
that this is where the two "fingers" of contamination 
caused by main sequence and background giant stars re- 
spectively, cross the cluster isochrone. Candidates from 
1.33 <V - Ic < 2.42 (15.5 <V < 18.0) suffer very little 
contamination, but note that this is largely because the 
addition oi B — V data removes the majority of the con- 
tamination in the VI CMD. For V - Ic > 2.42 {V > 18.0) 
there is growing incompleteness in the B band as well as a 
growing number of contaminants from the VI CMD alone, 
resulting in a drop in the discrimination in our member- 
ship selection. We expect that the level of contamina- 
tion is underestimated for the coolest bin (approximately 
V > 19.5) in Table ^, because the data below the clus- 
ter isochrone in the VI CMD is incomplete. The errors 
in the membership probabilities are estimated assuming 
Poisson errors in the numbers of candidate members, the 
numbers of candidate members rejected on the basis of 
their B ^ V and the numbers of predicted contaminants 
in the VI CMD. This latter error is an overestimate, be- 
cause the contaminant numbers arise from modelling pop- 
ulations several times larger. In any case, these errors are 
similar (as a fraction) to the simple Poissonian errors in 
the numbers of cluster candidates in each colour bin. In 
what follows we will simply assume that the contamina- 
tion fraction is known accurately, but will compare our 
results with what would have been obtained without any 
correction for contamination. 

It is of interest to compare the estimates above with in- 
formation appearing in the literature. Jeffries et al. (1998 ) 
used a similar BVI photometric selection technique to 
choose a cluster candidate list. These were subsequently 
followed up with spectroscopy. In that paper, 22 out of 31 
objects with 12.5 < V < 15.0 were confirmed as members 
based on their radial velocities. This is perfectly consistent 
with the estimates in Table |[ 

Hawley et al. (1999) did do an offset field (about 1 
degree from the cluster centre) covering 225 square ar- 
cminutes in V and / only. It is not clear (see Sect. [4.5D 
that this is far enough from the cluster centre to guarantee 
no cluster members and indeed Hawley et al. spectroscop- 
ically identified a few low mass NGC 2516 candidates in 
this field. However, we can use similar membership crite- 
ria to those in our survey to find how many bogus NGC 
2516 members we might expect in a 225 square arcminute 
area. This is then multiplied by 13.7 to match our survey 
area. The total number of contaminants expected in the 
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Table 4. Estimates of contamination and membership probability among our cluster candidates. Rows 1 and 2 list 
the observed numbers of cluster candidates in each colour interval, for single stars and binaries respectively. Rows 3 
and 4 list the predicted number of contaminants derived from our fitting technique in the VI CMD (see text). Rows 
5 and 6 list how many of these contaminants are rejected using the B — V colour and rows 7 and 8 give the final 
membership probabilities for candidates. 



V-I, 


0.53-0.70 


0.70-0.87 


0.87-1.05 


1.05-1.33 


1.33-1.68 


1.68-2.04 


2.04-2.42 


2.42-2.73 


2.73-3.02 


Observed 




















Single Candidates 


31 


61 


51 


115 


67 


66 


119 


146 


106 


Binary Candidates 


17 


41 


37 


62 


27 


23 


40 


40 


47 


Predicted 




















Single contaminants 


6 


38 


40 


63 


14 


19 


25 


63 


31 


Binary contaminants 


5 


18 


25 


55 


15 


11 


22 


39 


26 


Rejected 




















Single contaminants 


6 


22 


16 


32 


13 


13 


17 


26 


1 


Binary contaminants 


3 


8 


4 


27 


11 


11 


5 


19 


11 


Membership Prob. 




















Single candidates 


100 ±11 


74 ±13 


53 ±16 


73 ±9 


98 ±8 


91 ±9 


93 ±5 


75 ±7 


72 ±6 


Binary candidates 


88 ±17 


76 ±13 


43 ±17 


55 ±16 


93 ±19 


100 ± 20 


58 ±15 


50 ±20 


68 ±14 



VI CMD (to be compared with the sum of rows 3 and 4 of 
Table |) are 192 ± 51 for 0.87 < V ~ I^ < 1.68 (compared 
with 212 in Table|), and 274±61 for 1.68 < V-I^ < 2.73 
(compared with 179 in Table H). The former estimate is 
in good agreement with our own, but the latter is a little 
higher, perhaps indicating that there are indeed some low 
mass NGC 2516 members even 1 degree (^ 7pc) from the 
cluster centre. 



4. Luminosity and Mass Functions 

The CTIO catalogue of candidate members, together with 
the estimates of contamination in Table |4| allow us to 
make a preliminary investigation of the luminosity and 
mass functions (LF and MF) of NGC 2516 (that is the 
number of stars per V magnitude interval and the number 
of stars per logarithmic mass interval) . 



4.1. Luminosity Function 

The CTIO candidate membership catalogue is incomplete 
for V > 20, but is also incomplete for brighter stars be- 
cause of saturation in the CCD frames. As discussed in 
Sect. 2.4, there were 68 objects from the catalogue of 



bright stars in the NGC 2516 field compiled by Jeffries 
|et al. (1997 ), which are not present in the CTIO catalogue 
55 of these objects (which have V and B — V photome- 
try) reside within our CCD frames and 48 have photom- 
etry consistent with main sequence cluster membership. 
These should be i ncluded when calcu lating the LF and 
MF, and, because Jeffries ct al. (1997 ) estimate that this 
bright catalogue is complete to V^ = 13.5, we then have 
an almost complete catalogue of cluster members from 
5.8 < V < 20, covering a mass range of approximately 




Fig. 12. The V band luminosity function (LF) of NGC 
2516. The solid histogram shows the LF for candidate 
members identified from either B — V or V — Ic data 
(including brighter stars from Jeffries et al. 1997 ) and 



the dashed histogram represents the LF of those stars for 
which both B — V and V — Ic indicate membership. The 
two short lines at the base of the plot indicate the approxi- 
mate 90 percent completeness limits for these two samples 
at F = 20.0 and V = 18.2 respectively. 



0.3 < M < S.OMq (the upper limit being rather age, 
model and metallicity dependent). 
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13.5 <V< 15.5. An NGC 2516 LF corrected for con- 
tamination by non-members and plotted with M^ on the 
X-axis (assuming an apparent distance modulus of 8.3) ap- 
pears in Fig, [l^. We c ompare this with LFs from the field 
(|Reid fc Hawley 2000D and the Hyades (|Reid 1992|; |Reid fc 



Hawley 2000), which are normalised to match NGC 2516 



for 5 < My < 7. The correction for contamination in NGC 
2516 is applied by making each star worth < 1 in the lu- 
minosity histogram according to it's V — Ic colour and 
whether it is a photometric binary candidate, using the 
fraction of such stars which are likely to be members ac- 
cording to Table |4| For instance, a candidate member with 
V — Ic = 1.0 and which is flagged as a probable binary sys- 
tem will only contribute 0.43 to the LF bin corresponding 
to its V magnitude. All the very bright candidate mem- 
bers without V — Ic are assumed to be genuine members. 
The NGC 2516 LF in Fig. [l| exhibits both similarities 
and differences when compared with the LFs of the field, 
the Pleiades (which is almost iden t ical to the field LF 



Hambly, Hawkins, & Jameson 1993; Mcusinger, Schilbach 
fc Souchay 1996| ) and the Hyades. The LF is consistent 



with a monotonically rising curve at least as far as Afy — 
1 {V ~ 15.3 in NGC 2516). There are relatively more 
bright stars in NGC 2516 than the field because it is a very 
young cluster. There is then, despite the correction for 
contamination, a significant minimum at 7.7 < My < 8.7 
(16 <V < 17), that coincides with the "Widen dip" (see 



Wielen 1974; Upgren fc Armandroff 1981; BahcaU 1986) 



Fig. 13. The corrected absolute V magnitude luminosity 

f Unct ion (T,F - tjpp tPYt) nf KCC 9.^1 fi (gVinwn with prrnr 



that is also seen (at the level of 30-40 percent) in the LFs 
of the field, the Pleiades and the young a Persei cluster 
at My ~ 7.5 (grosser 1992|; |Mei.isinger et al. 1996|; [Belikov 



bars) ^nmparprl with T.Ft; fnr thp fiplH and tViP TTyarlpti 



The field and Hyades LFs are shown as histograms and are 
normalised to match the NGC 2516 LF at 5 < My < 7. 
The error bars on the field LF are about the same as in 
NGC 2516 at the same absolute magnitude, whereas the 
Hyades LF error bars are about a factor of two bigger (but 
with twice as many bins). The fie ld and Hyades LFs were 
taken from Reid fc Hawley (200C ) and references therein. 



et al. 199^; |Reid fc Hawley 2000| ) and in the Hyades at 
My = 8.5 ( Reid 1992| ). This dip is defined with reasonable 



Figure |l^ shows the LF of the candidate members 
without a correction for contamination. The plot shows 
the LF due to candidates selected using both B — V and 
V — I and also the complete candidate list (including the 
brighter stars discussed above and those stars classed as 
members for which either B — VovV — Ic are unavailable). 
We also mark on the plot the approximate faint magni- 
tude completeness limits for these two samples (discussed 
in Sect. 2.3). 

The rapid fall off at faint magnitudes coincides with 
our estimates of where the survey begins to be signif- 
icantly incomplete. There appears to be a minimum in 
the LF f or 15.5 < V < 17.0, w hich was also commented 
upon by Hawley et al. (1999| ). However, we might sus- 
pect that this feature merely results from the inclusion of 
many background contaminants at slightly brighter mag- 
nitudes (see Table ^ and Fig. ||) , resulting in the peak at 



clarity in NGC 2516, due to the large numbers of cluster 
members. There is some evidence that the overall rise in 
the LF towards fainter magnitudes then levels off for My > 
9 and certainly there is no sign of the very steep rise in the 
LF of a factor 2-3 between 9 < My < 10, that is seen in 
the Pleiades and nearby field stars. In this respect, NGC 
2516 is much more similar to the Hyades. 

We do not believe that incompleteness can be respon- 
sible for a lack of faint stars in NGC 2516. We have cal- 
culated that the LF is complete to at least My = 11.7, in 
the sense that stars are not missed because they were not 
detected. Perhaps then the LF of NGC 2516 is different 
to that of the Pleiades and the field, but there are other 
possibilities: First, it may be that that the contamination 
fraction is increasingly overestimated for y > 17 or un- 
derestimated for 14 < y < 17. We regard this as unlikely 
and in fact for the last column in Table 0, we believe the 
contamination fraction is probably underestimated (see 
Sect. [3.4| ). Second, we may have set our membership crite- 
ria too tightly and missed a significant fraction of fainter 
cluster members that lie either above or below our mem- 
bership bounds in the VI CMD. We already regard our 
membership criteria as quite generous and thus the fac- 
tor of approximately two increase in the LF that would 
be required between 10 < My < 12, seems unlikely to be 
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explained in this way. Third, it might be that mass seg- 
regation has been successful in removing lower mass stars 
from the central cluster regions and hence our survey. This 
has been put forward as an explanation of the differences 
between the Hyades and field LFs (Reid & Hawley 2000) 
and although NGC 2516 is much younger than the Hyades, 
we have only surveyed the central regions. Mass segrega- 
tion is discussed in more detail in Sect. 4.5. Lastly, it is 



possible that the mass function in NGC 2516 is similar 
to that in the Pleiades, even though the LF is not. This 
might be the case if the mass-luminosity relationship were 
markedly different at low masses because NGC 2516 has 
a low metallicity compare d wi th the Pleiades. This possi- 
bility is discussed in Sect. 4.2. 



a = +1,47+/-0.11_ 



a = -0.75+/-0,20 



E Sless et ol. Solar 



4.2. Mass Function 

The cluster MF can be explored in a preliminary way by 
calculating the mass of each candidate member using a 
mass-colour relationship derived from the isochrones we 
used to fit the cluster VI CMD. Unresolved binarity is of 
course a problem here. The net result of assuming that 
unresolved binaries are single stars, will be that the true 
MF will be slightly steeper than that derived (in the sense 
that an assumed single star of a given mass is actually two 
stars, one of which could have a much lower mass - see 
below). In principle, the numbers of low mass stars could 
be doubled by hiding them in binary systems with more 
luminous companions. 

The binarity in NGC 2516 and its effect on the de- 
rived MFs is discussed further in Sect. ^, but we note 
that in Sect. 3.2 we identified candidate members which 



-0.5 



are approximately equal mass binary systems based on 
their position above the cluster isochrone. Using the BV 
CMD we can also identify candidate binaries from among 



the bright s tars added to the CTIO sample from [Jeffries 



et al. (1997| ). Each identified binary system is counted as 
if it consisted of two identical components with the same 
colour. This is an approximation, because for objects with 
mass ratios less than unity, one component will be cooler. 
However, in the absence of better binarity information, 
a more complex approach is not warranted. Corrections 
for contamination by non-members are applied by using 
the fractional membership probabilities in Table |[ Each 
candidate will contribute a fraction to the MF bin cor- 
responding to its mass, according to the value oi V — Ic 
(or B — V a V — Ic is unavailable) and whether it is a 
candidate binary. 

The mass-colour relationships are derived from the 
same evolutionary models and empirical isochrones used 
to select the cluster candidates in the first place. The un- 
certain metallicity of NGC 2516 causes some difficulties 
here. Stars of a given colour have smaller masses at lower 
metallicities. We perform our calculatio ns using the so 



lar and half-solar metallicity models of giess et al. 2000 



and the solar metallicity model of D'Antona fc Mazzitelli 



(1997 ), which were discussed in Sect. \i.2\ Note that we 
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Fig. 14. The corrected mass functions for NGC 2516 (see 
text), derived from the solar metallicity modelofpiesE 
3t al. (200C|) (top), the solar metallicity model of D'Antona 
fe Mazzitelli (1997) (middle) and the half-solar metallicity 
model of ^iess et al. (2000| ) (bottom). The dashed lines 
indicate the power law fits discussed in Sect. 4.2 and the 



dash-dot lines indicate the approximate mass at which our 
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choose to use V — Ic, rather than B — V, as the primary 
indicator of T^g and hence mass, because it is relatively 
insensitive to metallicity (Leggett et al. 1996) and covers 
a wider range of masses. Only the hottest stars have their 
masses calculated from B — V, where i? — F is an almost 
metallicity-independent Tcff proxy (Castelli 1999|), so the 



and the results it yields in later sections are in reasonable 
agreement with the Siess et al. solar metallicity model. 

The power law fits to the higher mass range are in 
good agreement with the Salpeter value, agree well with 
the average value of a = +1.40 ±0.13 found for interme- 



fact thFct 



hg 



we iiave iguuicd tlic mutallicit^^ dcpunUeuce oi 
the B ~ V-Tci{ lelatiuuship in duiiviug the low-niutallicit3^ 
isochrone should be unimportant (see Sect. p^). 

The corrected MFs are shown in Fig. 113, expressed as 
the number of stars per unit logarithmic mass interval, as a 
function of logarithmic mass (in solar units) . The approx- 
imate lower mass limit at which our survey is complete 
(corresponding to F = 20) is about M = 0.3Mq for the 
solar metallicity models and M = 0.25Mq for the half- 
solar metallicity mo del. In these dia grams, the canonical 
stellar initial MP of ^alpeter (1955 ) would be a straight 
line of the form dN/d\ogM oc M"" with a = -hi. 35. 

Our MF for NGC 2516 is well defined between 3Mq 
(where any uncertainty in the age of NGC 2516 would 
start to be a factor and also where small number statistics 



diate mass stars in many open clusters by Phelps fc Jane 



(1993 ) and are close to the "universal" fie ld initial MF o f 



a = +1.S ± 0.3 for M > 0.5Mq found by [Kroupa (20011) . 
We also note that these quoted results neglect the effects 
of binarity and used a single relationship between V and 
mass to calculate the MF from the LF. This has the ef- 
fect of making the MF slightly less steep at high masses 
because binary systems are then treated as one star with 
a slightly higher mass. We have investigated what differ- 
ence this makes by simply treating our data in the same 
way. If we were to adopt a single relationship between V 
and mass, irrespective of binary status, then the slopes of 
our derived MFs would be smaller by about 0.2 and hence 
in even better agreement with the previously published 
values. 



b ecomj; iiupuiLaiiL) and 0.25-0.3Mq (wlieie iucunipleLe- 
ness Sets in). NuLe LliaL ueluw Lhe cuiiipleLenesb InniL we 
cannot simply say that the points are lower limits because 
there is also the uncertain level of contamination by cluster 



Determ inations of t he MF in the Pleiades ([Meusingei 



et al. 1996|) and M35 ( [Barrado y Navascues et al. 200lD 
clusters find MF slopes of a = -1-1.4 and a = -f 1.59 ±0.04 
for stars with M > 0.8Mq. Meusinger et al. then find that 



jLallici- 



non-m(|;mueis Lu euusiaei. i'ui eacli ui Llie Lwu mcL 
ties, tl leie is a clear Sdlpetei-like lise in the MF as the 
mass decreases, followed by a peak at Al ~ O.TA/q and a 
turnover towards lower masses. We have checked that the 
corrections for contamination by non-members and for bi- 
naries have very little effect on this overall shape. 



a c± in the Pleiades for 0.3-1.0Mq, whilst Hambly et al 
(1999D estimate a = -0.3 for M < 0.5Mq in the Pleiades 



Th(|; exact torm ot the Ml* might ot course be age 



and Barrado y Navascues et al. finds a = —0.2 for M35 
between 0.8 and O.2M0. The latter two slopes take no 
account of binarity and should probably be increased by 
~ 0.2 (see above) before comparison with our results. At 
very low masses it is likely that the MF falls again. Bouviei 



model and metahicity-dependent. barrado y Navascues 
et al. |200l[) investigated the MF of M35, a rich north- 
ern hemisphere cluster of similar age to the Pleiades and 
NGC 2516, using a variety of ages and (solar metallic- 
ity) evolutionary models. Their results show that above 
0.2Mq, derived MFs are virtually identical in all cases. 

We have parameterised our derived MFs in terms of 
two power law fits in the range 0.3-0.7Mq (or 0.25- 
0.7Mq for the half-solar metallicity model) and 0.7- 
3.OM0. For the higher mass range we find a = ±1.47 ± 
0.11 and a = ±1.67 ± 0.11 for the solar and half- 
solar metallicity models of ^iess et al. (2000| ), and a = 
±1.58 ± 0.10 for the solar metallicity model of pD'Antona 
& Maizitelli (1997). For the lower mass range we find 
a = -0.75 ± 0.20 and a = -0.49 ± 0.13 for the so- 
lar and half-solar metallicity models of Siess et al., and 
a = —1.00 ± 0.18 for the solar metallicity model of 
D'Antona & Mazzitelli. The systematic difference between 
the models is caused by the change in the colour-mass re- 
lationship. At lower metallicities, stars of the same colour 
have lower masses, increasing the value of a, but there 
is also some modest model-dependence. For simplicity of 
discussion in what follows we shall only use the two Siess 
et al. models. We have also done all our calculations for 



et al. (1998D and [Moraux, Bouvier, fc Stauffer (200l| ) find 
a ~ —0.4 and —0.5 across the brown dwarf boundary at 
^ 0.05-0.2Mq in the Pleiades and Barrado y Navascues 
et al. estimate a more extreme slope of —1.8 below 0.2Mq 
in M35. This behaviour is mirrored in the field MF, where 
a Salpeter-like slope is found above 0.6Mq, a relatively 
flat MF with an a of -0.1 to ±0.3 down to - O.IMq, 
and then a decline into the brown dwarf regime with a 
between -0.5 and - 1.0 (|Gould, Bahcall, fc Flynn 199*7 ; 



Rcid et al. 1999; Chabrier fc Baraffe 200C; Kroupa 2001) 



If NGC 2516 has a solar metallicity, then its MF drops 
much more sharply towards lower masses than occurs in 
either the Pleiades or field for 0.3 < M < O.7M0. This 
steep slope is of course directly related to the deficit of 
low luminosity stars in the NGC 2516 LF with respect 
to the Pleiades and the field which we remarked upon in 
Sect. 4.1. Any of the explanations offered there (such as 



mass-segregation) might make the drop in the MF below 
0.7Mq less extreme. A Pleiades-like MF could still be 
recovered for the solar metallicity scenario if the number of 
NGC 2516 members between 0.3 and O.5M0 were roughly 
doubled with respect to the more massive stars. If NGC 
2516 has a half-solar metallicity then the discrepancy with 
the Pleiades and field MFs is reduced but still present. 



the solar metallicity [D'Antona fc Mazzitelli (1997| ) model 
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Fig. 15. Binary fraction in NGC 25 16 for systems with 
approximately 0.6 < g < 1 (see Sect. O) as a function of 
V-L. 



4.3. Binarity 

The frequency of binary systems and the distribution of 
their mass ratios are important constraints on star forma- 
tion theories. So far in this work we have made a partial 
correction for binarity in determining the cluster MF, by 
identifying cluster candidates with near equal mass ratios 
{q = M2/M1 ~ 1). Unresolved binarity will also influence 
the deduced total mass of the cluster (see Sect. i.4). We 
have also not considered the possibility of triple systems 
in this paper, which could place an object up to 1.2 mag. 
above the cluster isochrone (for equal mass components). 
However, the numbers of missed triple systems should be 



allow these isochrones to move up or down by ±0.1 mag- 
nitude (see Sect. [3.2| ), the mean binary fraction changes 
by ±5 percent. As we discussed in Sect. 3.2, the shape of 
the half-solar metallicity isochrone is slightly different for 
V — Ic > 2.0, which would result in a differential increase 
in the binary fraction of cool stars with respect to the hot 
stars of about 5 percent. This would be enough to remove 
the possible trend just discussed, so it is premature to 
claim to have seen any dependence of binary fraction on 
mass. 

The binary fraction estimated in this way is a lower 
limit. The mass ratios to which our photometric identifi- 
cation technique is sensitive, ca n be estimated using the 
results and equations found in Kahler (1999 ). We note 
that the q sensitivity limit is not quite independent of 
colour for stars 0.3 mag above the single star locus, be- 
cause the gradient of the V vs V ~ Ic cluster isochrone is 



not constant. Using equation 5 in Kahler (199£ ) and the 
gradients determined from the cluster isochrone we find 
that we are sensitive to 0.59 < q < 1.0 for V — Ic < 0.6, 
0.55 < g < 1.0 for 0.5 < V - Ic < 1.1, 0.62 < g < 1.0 for 
l.KV ~ Ic < 2.5 and 0.60 < q < l.O for V - Ic > 2.5. 
Given an approximately fiat distribution of q or one that 
rises towards lower q values, any small variations in appar- 
ent binary fraction seen in Fig. uB could also be partially 
explained by this varying sensitivity. 

In field stars, Duquennoy & Mayor (1991) found that 
the mass ratio distribution increased towards lower mass 
ratios, perhaps peaking at g = 0.3, and that the binary 
fraction of solar-type stars was about 70 percent for all 
mass ratios, and 21 percent for 0.6 < q < 1.0. Work on 
the Pleiades s olar-type stars is in broad agreement with 
these results ( Mermilliod et al. 1992 ; Bouvier, Rigaut, fc 



quite rfire - less than 2 percent of systems in the Pleiades Nadeau 1997|) . There is some evidence that the binary 



( [Mermilliod et al. 19921 ). 

Figure |15| shows the fraction of photometrically iden- 
tified binaries as a function oiV — Ic colour in NGC 2516, 



where nembership was determined from the solar metal- 



licity isochrone of Siess et al. (200C ). The binary frac- 
tion is defined as the number of probable binary systems 
divided by the total number of systems. The fractions 
are calculated using our membership catalogue, the bi- 
narity fiags and the membership probabilities in Table ^. 



fraction among lower mass field objects is smaller (30-40 
percent), but that they are more inclined to be found in 
q 2± 1 systems ( Fisher fc Marcy 1992 ; Reid fc Hawley 
2000|) . Using a simil ar photometric selecti on technique to 



the one used here, ptauffer et al. (1984]) found that 26 
percent of Pleiades low mass stars (0.7 < y — / < 2.1) 
were more than 0.3 mag. above a single star isochrone and 
that this was similar to the 22 percent photometrically 
determined binary fraction in hotter stars found by Bettis 



The weighted mean binary fraction is 26 ± 2 percent with 
marginal evidence (at the 90 percent confidence level) for a 
decrease in the binary fraction towards lower masses. This 
may be more significant, as we expect the lowest mass bin 
in Fig. na to have an overestimated binary fraction, be- 
cause incompleteness sets in and will bias against single 
star detection. It is also the case that low mass M-dwarfs 
are more prone to flaring, which might result in some frac- 
tion of the lowest mass single stars being misclassified as 



(1975) 



The main result of this subsection is that the binary 
fraction in NGC 2516 for stars with q in the approxi- 
mate range 0.6-1 is similar to that of the Pleiades or the 
field, when the binaries are identified in the same way. If 
the distribution of q were that proposed for the field by 
Duquennoy & Mayor (1991), down to q = 0, then the total 



binaries (see Stauffer et al. 1984). The confidence level for 
a significant downward trend increases to 98 percent if 
this last point is excluded. On the other hand, the error 
in the binary fraction should be increased to take into ac- 
count the somewhat arbitrary nature of the isochrones we 
have used to define the cluster single star sequence. If we 



binary fraction in NGC 2516 would be about 85 percent. 
If the distribution of q were flat, then the total binary 
fraction would be 65 percent. 

A subsidiary consideration is estimating what effect 
unresolved binarity might have on the derived mass func- 
tions. Our identification oi q > 0.6 binaries, where the 
companion star can significantly change the colours and 
magnitudes of the system as a whole, results in MFs which 
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have slopes about 0.2 larger than those where this is ig- 
nored (see Sect. 4.2). These MFs have the merit of beoing 
obtained from observables which are in principle compara- 
ble with other estimates derived fr'om photometry data in 
other clusters, and do not rely on knowledge of the mass 
ratio distribution. However, even for those systems with 
q < 0.6 there can be a significant number of "hidden" 
lower mass stars, resulting in a MF for all stars which 



is even steeper still. Bagar & Richtler (1991) investigated 
this effect in 2-14Mq stars in LMC clusters. They found, 
for true MFs with a = -1-1.5, a binary fraction between 0.5 
and 1, and binary companions drawn randomly from the 
same MF, that the observed MF (estimated from a mass- 
magnitude relation) had a slope that was 0.3-0.4 smaller. 
The effect was even larger for smaller values of a. A true 
a of -1-0.5 might appear as a = — 0.4 for a bi nary fraction 
of 1.0. Similar simulations by Kroupa (2001 ) yield differ- 
ences between the system and single star MFs largely in 
agreement with these results. 

We point out to the reader that by using a colour-mass 
relation and identifying and dealing with q > 0.6 binaries, 
we have partially alleviated this problem. To try an gauge 
by how much the true stellar MF slopes might be further 
increased over the quoted MF slopes in this paper, we 
randomly added binary companions to a fraction of the 
"single" stars. Using a total binary fraction of 65 percent 
and a flat q distribution we find that the true single-star 
MF a is increased by a further 0.05 for 0.7 < M < 3.0Mq 
and 0.3 for 0.3 < M < O.TMq (using the piess et al. 
2000rsplar metallicity model s). Using the q distribution 
proposed for field binaries by Duquennoy fc Mayor (1991 ) 
and assuming a total binary fraction of 85 percent, results 
in a increasing by 0.1 and 0.4 in these two mass ranges. 



4.4. The cluster mass 

We can integrate our corrected MF (from the solar metal- 
licity isochrone) to yield an estimate of the total mass 
of the cluster (down to about 0.3Mq) of IIOSMq. More 
than half of this mass is contained in stars with 0.6-2Mq 
with decreasing contributions at lower and higher masses. 
The equivalent calculation for the MF derived from the 
half-solar metallicity isochrone is 945Mq (complete to 
O.25M0). 

Unresolved binarity will increase the derived cluster 
mass. We have made a partial correction for this in our 
work so far, by identifying near equal mass-ratio binary 
systems. If we had not done so, our deduced cluster mass 
would have been only SSOMq (solar metallicity). To es- 
timate the maximum likely contribution that could arise 
from binaries with q < 0.6, we can use the q distribution 



proposed for field binaries by Duquennoy & Mayor (1991) 
and assume that the total binary fraction is 85 percent 
and independent of primary mass. Integrating this distri- 
bution, we find that the cluster would be 35 percent more 
massive than if all the stars were single, and thus the to- 
tal cluster mass would be 1190Mq (solar metallicity), for 



primaries with M > 0.3Mq. A correction could also be 
applied for stars less massive than this. Integrating an ex- 
trapolated MF derived from the 3iess et al. (200C) solar 
metallicity isochrone (see Sect. 4.2) from say 0.0-0.3Mq, 
yields only another 73Mq. The corresponding additional 
mass for the half-solar metallicity model is 72Mq (0.0- 
0.25Mq). Thus irrespective of the cluster metallicity, the 
contribution of low-mass stars and brown dwarfs to the 
total cluster mass inside our surveyed area is likely to be 
less than 10 percent unless there were a sharp upturn in 
the MF below O.3M0. 

The cluster mass is similar, but a little higher, than 
that found for the Pleiades in the same mass range. 
Meusinger et al. (1996) quote a figure of 800Mq for 
M > 0.3M(i^ and [Pinficld et al. (1998] ) derive 735M0 
for all masses down to the substellar limit. However, we 
have to be careful to compare like with like. The quoted 
Pleiades results are extrapolations for all stars out to the 
cluster tidal radius. 

Pinfield et al. (1998) give an approximate expression 
for the tidal radius of a cluster in a circular Galactic orbit, 
close to the Sun as 



n 



GAL 



2{A-B) 



1/3 



(4) 



where Mc is the cluster mass and A and B are the Oort 
constants in the solar neighbourhood. Using the value of 
A — B from Kerr & Lynden-Bell (1986), this reduces to 



1 /3 

rt — 1. 4:6 Mc , with M^ in solar masses and rt in parsecs. 
If we take IOOOMq as the minimum mass for NGC 2516, 
rt > 14.6 pc. 

Using a distance modulus of 7.9, our survey of the clus- 
ter covers a square area, 6.2 pc on a side. Thus we expect 
significant numbers of cluster members outside the area 
covered by our survey. Pinfield et al. (1998) find that the 
Pleiades has rt — 13.1pc and that although the high mass 
stars are highly concentrated within a few pc of the clus- 
ter centre, there are significant numbers of low mass stars 
at much greater distances, such that only half the cluster 
mass was contained within 3.66 pc. If NGC 2516 is analo- 
gous to the Pleiades in terms of the spatial distribution of 
its members and the amount of mass segregation present 



(see Sect. 4.5), then the total cluster mass might be sig- 
nificantly greater than just what we have observed in our 
limited survey. 

4.5. Mass segregation 

Dynamical evolution and mass segregation can have a sig- 
nificant effect on the shape of the present day LF and 
MF of NGC 2516. Equipartition leads to a growing core 
radius with decreasing mass and the preferential eva pora- 
tion of low-mass cluster m e mbers (see for example pe la 
Fuente 1995; Kroupa 1995; ie la Fuente Marcos fc de la 



Fuente Marcos 2000| ; [Kroupa, Aarseth, fc Hurley 200l|) . 
The evidence for dynamical effects on the MFs of open 
[clusters older than a few hundred Myr is strong (e.g. Reid 
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text). 

ulus of 7.9. The horizontal dashed lines indicate the level 
of (assumed uniform) contamination expected in each of 
these samples. The solid lines are King profile fits for as- 
sumed tidal radii of 14.6 pc (see Table ||). 



1992; [Montgomery, Marschall, fc Janes 1993| ; parajcdini 
et al. 1 999 ). The segregation effects are less obvious, but 



s till seen, in younger clusters like the Pleiades and M3 5 



(Pinfield et al. 1998; Barrado y Navascues et al. 2001) 



the lowest mass stars, perhaps because these clusters have 
ages similar to their dynamical relaxation timescales. 

As a step towards understanding how much mass seg- 
regation might have taken place in NGC 2516, we exam- 
ine the radial surface density profiles of candidate clus- 
ter members split up into bins according to their V — Ic 
colours and hence masses. We exclude binary candidates 
at this stage, because their colours are not reliable mass 
indicators - they could be up to twice as massive. We 
have to acknowledge though that we must have included 
some binary stars with q < 0.6, which will have the effect 
of blurring any distinction between different mass bins, 
perhaps lessening the effects of mass segregation. In or- 
der to have sufficient stars to draw statistically signifi- 
cant conclusions we divide the stars into just 4 mass bins 
(see Table g|). The division is made preserving the colour 
boundaries in Table ^ so that estimates of field star con- 
tamination can still be made. We do not include the lowest 
mass interval from Table ^ because of possible uncertain- 
ties in the exact level of contamination (see Sect. 3.4). 
Note that the values for the boundaries of the mass bins 
are derived from the solar metallicity 3iess et al. (2000|) 
isochrone. The corresponding mass boundaries in Table l5| 
for the half-solar metallicity would be: 1.38Mq, 0.79Mq, 
O.58M0 and O.27M0. 

The centre of mass for each bin is determined by min- 
imising the projected moment of the stars about a point. 
The centres drift south-west by about 2.5 arcminutes be- 
tween the highest and lowest mass bins. This small drift 
has no effect on our results, which would be almost iden- 
tical if we had fixed the cluster centre at a single position. 
The surface density of stars as a function of radius (assum- 
ing a cluster distance modulus of 7.9) is then calculated 
for each mass bin, taking account of the geometry of the 
surveyed area. 

The results are shown in Fig. |lj. The horizontal dashed 
line in each plot indicates the (assumed uniform) density 
of contaminating objects. It now becomes clear that the 
probability that a candidate member is actually a con- 
taminating object shrinks drastically if we only consider 
objects within ^ 15 arcminutes {^ 1.7pc) of the cluster 
centre. Conversely, if we consider objects in the outer part 
of our survey, the probability of membership is lower than 
the average values in Table |j. 

To get some parameterisation of the mass segregation 
we have fitted empirical profiles of the form (see King 



v-a.vis valiips assnmps a rliistpr distancn mnrl- Wb2\) 



P = Pa 



v/1 + {r/rcy VI + i^t/rcy 



(5) 



with a hint that equipartition has not been achieved in 



where po is a normalisation constant, Vc is the "core ra- 
dius", where the surface density falls to half its central 
value, and rt is the tidal radius where the cluster is trun- 
cated. The King model assumes that there is no contam- 
ination in the sample, so we subtract the levels indicated 
by the dashed lines in Fig. |l^ before fitting. Since our data 
are confined to regions well within r^ , the fits are relatively 
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Table 5. King profile fitting results (assuming a cluster distance modulus of 7.9). Columns 1 and 2 show the mass 
subset considered (for the solar metallicity isochrone - see text), col. 3 is the assumed tidal radius, col. 4 the fitted 
core radius, col. 5 the fitted normalisation constant, col. 6 the expected surface density of contaminating field stars 
and col. 7 the x^ of the fit (with 10 degrees of freedom). For each mass subset, cols. 8 and 9 list the total numbers of 
cluster members in our survey and the total number of cluster members out to the tidal radius predicted by the fitted 
model. 



V-/c 

subset 
<0.53 



Mass 

(Mq) 

> 1.48 



0.53-1.05 0.87-1.48 



1.05 - 1.68 0.65 - 0.87 



2.73 0.35 - 0.65 



Centre 
(J2000.0) 



RA = 
Dec = 
RA = 
Dec = 
RA = 
Dec = 
RA = 
Dec = 



7 58 03.2 
-60 45 29 
7 58 05.6 
-60 46 48 
7 57 56.0 
-60 47 00 
7 57 48.0 
-60 47 00 



rt 



14.6 
18.4 
14.6 
18.4 
14.6 
18.4 
14.6 
18.4 



(pc) 



0.841":- 



+0.16 

0.14 
+0.23 

0.17 
+0.21 

0.16 

Ql+0.56 

^^-0.40 

+0.53 

0.37 
+0.99 

0.64 
+0.87 

0.56 



0.81 

0.90 

0.86 

1 

1.80 

3.52 

3.26 



po 



Pcont 



(stars/arcmin ) 



0.223 
0.218 
0.228 
0.225 
0.141 
0.132 
0.279 
0.238 



0.0 
0.0 
0.013 
0.013 
0.010 
0.010 
0.017 
0.017 



X 



7.7 
7.8 
10.9 
11.0 
11.9 
11.9 
13.4 
13.1 



-^ 'SUr' 



92 
92 
103 
103 
150 
150 
280 
280 



iVtid 



113 
121 
126 
135 
196 
216 
637 
725 



insensitive to this parameter. To test this we performed 
two sets of fits, fixing r^ at either 14.6 pc (corresponding 
to a minimum cluster mass of IOOOMq - see equation 0) 
or 18.4 pc for a cluster mass of 2000Mq, which acknowl- 
edges that a significant fraction of the cluster mass may 
lie outside the CTIO survey. 

The fits are all reasonable, with x^ values of between 
7.7 and 13.4, with 10 degrees of freedom in each case. 
The fit results, including 68 percent confidence intervals 
in Tc, and the assumed surface density of contaminants, 
are reported in Table o. 

The King profile modelling suggests that the core ra- 
dius undergoes a dramatic enlargement as the mass de- 
creases from ~ I.OA/q to ^ 0.5A'/q and that this is in- 
dependent of the assumed value of rt. It could be argued 
that this might be explained in terms of an underesti- 
mation of the contamination level in the two lower mass 
bins. Such an underestimation would lead to a "flatter" 
distribution of surface density that would then be fitted 
with a large value of re- To test this, we re-fitted the 
two lowest mass bins with a King profile plus a variable 
uniform surface density, which simulates an extra unac- 
counted for level of contamination. The core radius was 
fixed at 0.9 pc, corresponding to the value of re for stars 
with mass 0.86 < M < 1.48Mq. We found that the x^ 
values were similar with this model (x^ values of 13.6 and 
1 1. 5 respectively for the same numbers of degrees of free- 
dom) , but the levels of contamination implied by the fits 
(i.e. the extra constant surface density ter ms) were factors 
of 3-5 higher than we estimated in Sect. 3.4, This would 
mean that our estimates of the contamination levels were 
wildly in error and that only about 35 percent of the clus- 
ter candidates with 1.05 < V — Ic < 2.73 were members 
in contrast to our present estimate of 84 percent. We do 
not believe that errors of this extent are possible. 

We conclude that the evidence for mass segregation 
is strong when comparing stars above and below about 
O.SMq. This coincides with the break in slope of the MF 
and lead us to suspect that the downturns in the LF and 



MF of NGC 2516 at low masses may be explained by mass 



segregation (see Sects. 4.1 and 4.2). We can test this by 
extrapolating and integrating the surface density profiles 
beyond the extent of the current survey. King (196^ gives 
an expression for the integral of equation |5| out to the tidal 
radius, which we use to calculate the total number of clus- 
ter stars in each mass bin, A'tidai, in Table ^. This num- 
ber can be compared with the actual number of cluster 
members (corrected for contamination), A'survcy, to esti- 
mate what fraction of cluster stars lie outside our survey. 
Naturally, such an extrapolation is reliant on a knowledge 
of the surface density profile outside the area we have sur- 
veyed! 

The fraction of cluster stars inside our survey de- 
creases from 81 percent for M > 1.48AfQ to 44 percent 
for 0.37 < M < O.64M0, if rt = 14.8 pc. These fractions 
decrease to 76 percent and 39 percent for the larger value 
of the cluster mass and rt . 

There are two major implications of this result. (1) The 
total cluster mass is probably greater than we first esti- 
mated. Most of the cluster mass is concentrated in stars 
between 0.6-2Mq so it is likely that the cluster is more 
massive by about a factor of 1.3. (2) The decrease in the 
MF towards lower masses will be flattened off. Earlier we 
stated that an increase of roughly a factor of two in the 
numbers of stars below 0.5Mq (with respect to the num- 
bers of more massive stars) could bring the LF of NGC 
2516 roughly into agreement with that of the Pleiades. 
The results in Table o show that within our survey area 
the ratio of stars with 0.37 < M < 0.64Mq to those with 
M > 1.48Mq is 3.0, but that this ratio for all cluster stars 
is probably as large as 5.6-6.0; approximately the increase 
required. Taking the results from the lowest two mass bins, 
and assuming that the correction to the MF will vary lin- 
early (with log mass) from about 1.3 at0.7MQ to 2.3 at 
0.3Mq, the fitted values of a (see Sect. 12) would increase 
from -0.75±0.20 to -0.08±0.19 and from -0.49±0.13 to 
-|- 0.13ib0.13 for the solar and half-solar metallicity models 
of Siess et al. (2000 ) respectively. 
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It is difficult to provide precise corrections to the MF 
and total cluster mass because our analysis is dependent 
to a great extent on the assumed form of the surface den- 
sity distribution outside of the region surveyed and to a 
lesser extent on the assumed non-member contamination 
fraction in each bin. For example, if the level of contami- 
nation for low-mass stars is higher than we have assumed 
then: the numbers of low mass stars inside our survey 
will be smaller; the fitted core radius will decrease and 
the fraction of cluster members inside the survey will in- 
crease ; the cluster mass inside the survey will be smaller; 



the Mi p will drop more steeply and the total cluster mass 
will be smaller. To give an idea of the possible size of these 
effects we ran through our analysis again, but systemat- 
ically increased the number of contaminants predicted in 
rows 3 and 4 of Table Q, for the columns with V — Ic> 1.68 
(M < 0.64Mq), by twice their Poisson errors. The results 
from this were: that the MF slope below 0.6Mq changed 
from a = —0.75 to a = —0.92; the cluster mass inside the 
survey decreased slightly from 1105Mq to IOSOMq; the 
core radius in the lowest mass bin decreased from 3.5 pc to 
2.8 pc; and the fraction of stars with 0.37 < M < 0.64Mq 
inside our survey increased from 44 to 52 percent. Changes 
of this order do not affect any of our major conclusions. 



5. Discussion 

NGC 2516 has been referred to as the "southern Pleiades" 
on the basis of thei r similar age s and apparent richness 
of high mass stars ( Eggen 1983 ). The central questions 
which we can attempt to answer on the basis of the results 
presented here are: (1) how similar are NGC 2516 and 
the Pleiades, in terms of their LFs, MFs, total mass, mass 
segregation and binarity, especially among low-mass stars? 
(2) What are the prospects for finding even lower mass 
objects and brown dwarfs in NGC 2516? 



The Pleiades survey of Meusinger et al. (1996) cov- 
ered the central 16.5 square degrees (about 9.3 x 9.3 pc) 
and was complete to M„ — 12. The LF shows a rela- 
tively sharp increase of a factor ~ 2 at My = 10 and 
the MF is quite flat, with a ~ (see Sect. 4.2), between 
0.3 < M < I.OMq. The LF and MF were simply corrected 
by a linear factor for stars of all masses and luminosities 
to account for Pleiades stars which lay outside the survey 
area. Thus no account of mass segregation was taken other 
than a crude correction to the cluster mass, which resulted 
in Mpiciadcs = 800Mq down to stars of 0.3Mq. Hambly 
et al. (1999| ) derived the MF from a wider 36 square de- 
gree survey (13.7 x 13.7 pc) and obtained a ~ —0.3 below 
O.5M0. 

Using survey ma terial complete to radii of about 6 pc, 
Pinfield et al. (199q ) show that core radii in the Pleiades 
increased from about 1 pc for stars of 3-5A/q solar masses 
to 3pc at 0.3Mq. As a result we can say that Meusinger 
et al. 's work must have been affected by mass segregation, 
with roughly one third of the low-mass members lying 
outside their survey, but essentially all the Pleiades stars 
with M > 1.2Mq included. Hambly et al.'s work fares 
better, with perhaps only ^ 20 percent of stars with M < 
0.5Mq missing. In both cases then, the true value of a is 
likely to be a little higher than derived from these spatially 
limited surveys, although the effect should not be nearly 
as large as calculated for NGC 2516 in Sect. 4.5. Pinfield 
et al. estimate that the total mass of the Pleiades, after 
accounting for this segregation is 735 Mq. 

The core radius values in the Pleiades are comparable 
with the NGC 2516 values that we determined in Table ra, 
but it does seem that the growth of re with decreasing 
mass is more rapid in NGC 2516. A simple power-law fit 
indicates that re oc M~^ with /3 ^ 0.9 ± 0.2, rather tha n 
the 0.5 found for the Pleiades by [Pinfield et al. (1998|) . 
In a virialised system, w e might expect Tc oc M~^'^ 
(see Pinfield et al. 199g for an analytic argument and 
ppitzer fc ShuU 1975 for numerical simulations). This is 



5.1. NGC 2516 and the Pleiades 

Of key importance in any comparison with the Pleiades is 



perhaps a hint that at least some of the mass segrega- 
tion in NGC 2516 is pr imordial, as has been supposed for 
some younger clusters ( ^agar et al. 1988 ; Bonncll fc Davie; 



to mak |e sure we arc comparing similar surveys, in terms of 1998]) , with the high-mass stars being initially more cen 



their area and completeness. We have seen that our NGC 
2516 survey is complete to T^ = 20, (M„ ~ 11.7) and cov- 
ers the central 6.2 x 6.2 pc. The Pleiades is approximately 
3 times closer and thus surveys of the Pleiades must cover 
9 times the area of our NGC 2516 survey to be equivalent. 
Dachs & Kabus (1989) presented the LF for bright 
stars {My < 5) in NGC 2516, from a 1 degree diameter 
central region. They showed that this LF was similar in 
shape to that of the Pleiades with approximately a fac- 
tor of two more early-type stars in NGC 2516. There was 
marginal evidence for mass segregation in their sample, 
but few stars at this luminosity were estimated to lie out- 
side their survey. On this basis they claimed that NGC 
2516 was about twice as massive as the Pleiades, implic- 
itly assuming that the two clusters have similarly shaped 
MFs at all masses. 



trally concentrated than the low-mass stars. This would 
be a mildly surprising result, because dynamical mass 
segregation should have removed the signature of initial 
conditions on the cluster relaxation timescale - which is 
about 100 Myr for clusters of the size of NGC 2516 and 
the Pleiades. 

Both the total mass and the MF of NGC 2516 are 
dependent on what is assumed for the metallicity of the 
cluster and also to a lesser extent on which stellar evolu- 
tion models are used. If the cluster has a solar metallic- 
ity then the mass inside our survey area (for stars with 
M > 0.3Mq) is 1100-1200Mq, depending on how un- 
resolved binarity is treated. The equivalent figure for a 
half-solar metallicity model is 95O-1O5OM0 (for stars with 
M > 0.25Mq). If the cluster does have a solar metal- 
licity then the MF drops sharply below 0.7Mq, with 
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a ~ —0.75. A half-solar metallicity yields a shallower 
slope of a ~ —0.49. We have established that most high 
mass (> 1.5Mq) stars are included in our survey, but that 
more than half of the low mass (< 0.6Mq) stars of the 
cluster may lie outside this region - depending on the ex- 
act form of the density distribution beyond our surveyed 
area. Correcting the total cluster mass for this segrega- 
tion leads to an estimate of about 1240-1560AfQ - twice 
the mass of the Pleiades, in agreement with the earlier 
prediction made by Dachs & Kabus (1989). This implies 



that the shapes of the MFs of NGC 2516 must be rea- 
sonably similar (at least over the mass range which con- 
tributes mos t to the total cluster mass). As we discussed 
in Sects. 4.2 and l4.5|, a factor of two increase in the num- 



bers of low mass stars relative to high mass stars could 
increase a at low masses and bring the NGC 2516 MF 
into into agreement with the flat (a ~ 0) Pleiades MF 
between 0.3 < M < I.OMq. Mass segregation appears to 
provide just this correction, increasing a by about 0.7. 

Recent simulations of an evolving star cluster by 
Kroupa et al. (2001 ) show both the total MF and the 



MF for stars within 2 pc of the centre of a cluster similar 
in size to the Pleiades, and at an age of 100 Myr. Their 
Figures 14 and 15, demonstrate that mass segregation can 
indeed produce a downturn in the MF of the central re- 
gions below 0.7AfQ, even if the whole-cluster MF is flat. 
A thorough survey for low-mass stars outside the area dis- 
cussed in this paper will be vital to constrain the surface 
density of members, the whole-cluster mass function and 
the total cluster mass. 

On the basis of a comparison of high mass ratio {q > 
0.6), unresolved binary systems, the Pleiades and NGC 
2516 have a similar binary fraction of 26 ±5 percent. It re- 
mains to be seen whether this similarity persists for lower 
mass ratio binary systems. 

In summary, subject to a wider survey confirming the 
presence of an extended population of low-mass objects in 
NGC 2516 (see the offset field of Hawley et al. 1999 for 
some preliminary evidence), we believe that NGC 2516 
and the Pleiades are similar in most respects, but that 
NGC 2516 has twice the mass and numbers of stars. 



5.2. Lower mass stars and brown dwarfs in NGC 2516? 

Any apparent difference between the Pleiades and NGC 
2516 MFs at low masses could be very significant for those 
searching for even lower mass stars and BDs in NGC 2516. 
If we take our observed NGC 2516 solar metallicity MF 
below 0.7Mq, with a — —0.75 (see Sect. |4^ ), and simply 
extrapolate to lower masses, we find that there would only 
be about 100 brown dwarfs (with 0.030 < M < 0.075Mq) 
in the area we have surveyed. The half-solar metallicity 
MF, with a — —0.49, predicts around 220. The total num- 
ber of brown dwarfs in the cluster could be much larger if 
mass segregation continues to lower masses. At the very 
least these numbers should be doubled to account for the 
segregation that is already apparent in ^ 0.5Mq stars. 



If instead we assume a Pleiades-like MF that is flat 
between 0.7Mq and 0.2Mq and then decreases with a = 
—0.5 below this, we calculate that there would be 360- 
440 brown dwarfs in the solar and half-solar metallicity 
cases respectively. Even with this MF, the total mass in 
the form of brown dwarfs would be less than 3 percent of 
the cluster mass and the extra contribution in the form 
of stars with M < O.3Af0 would be about 15 percent. A 
search for brown dwarfs and low mass stars in the outer 
regions of NGC 2516 will tell us much about the whole- 
cluster mass function and the cluster dynamics. 

6. Summary 

In this paper we have presented a large, accurate and uni- 
form BVIc survey of stars covering 0.86 square degrees 
(~ 6.2 X 6.2 pc) of NGC 2516, and which is almost com- 
plete to y = 20, M ~ O.3M0. Thanks to the relatively 
low contamination of the cluster main sequence by fore- 
ground and background objects we have been able to se- 
lect a sample of candidate cluster members. NGC 2516 
probably has a lower metallicity than the similarly aged 
Pleiades and our list of candidate members list will be 
an important source of optically selected low-mass targets 
for further investigations of magnetic activity, elemental 
abundances, lithium depletion and rotation rates in con- 
vective stars. 

We have made a preliminary investigation of the lu- 
minosity function, mass function, binarity, mass segrega- 
tion and total mass of NGC 2516, and compared our re- 
sults with the Pleiades, other young clusters and the field. 
Because the metallicity of NGC 2516 may be as low as 
half the solar value, we have performed two sets of calcu- 
lations, one for a mass fraction of heavy elements, z — 0.02 
and another with z = 0.01. Some of our results are some- 
what dependent on the assumed metallicity and others 
depend on extrapolating to lower masses or larger spatial 
areas. We can summarize the results of our investigation 
as follows: 

1. The luminosity function of the central regions of NGC 
2516 are consistent with the Pleiades and field popu- 
lations for Aft, < 8. For fainter stars, the luminosity 
function of NGC 2516 flattens and contains a factor of 
two fewer stars than are seen in the Pleiades and field. 

2. The derived mass function for NGC 2516 is in agree- 
ment with the Pleiades and field for stars with 0.7 < 
M < 3.0A^Q, with some metallicity and evolution- 
ary model dependence. We find a mass function, 
dN/d\ogM oc Af-°, with a = -hl.47 ± 0.11 and 
a — +1.67±0.11 for the solar and half-solar metallicity 
models of piess et al. (2000] ) and a = -Hi. 58 ± 0.10 for 



the so lar metallicity models of D'Antona fc Mazzitell: 
(1997 ) . These values are close to the canonical ^alpetei 
(1955 ) value of a = -1-1.35. At lower masses (0.3 < 
M < 0.7Mq) we find a = -0.75 ± 0.20 and a = 
—0.49 ± 0.13 for the solar and half-solar metallicity 
Siess et al models, and a — —1.00 ± 0.18 for the solar 
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metallicty D'Antona & Mazzitelli models. These are 
more steeply declining functions of mass than seen in 
either the Pleiades or field. 

3. Mass segregation is clearly present in NGC 2516 when 
we compare the radial distributions of stars above and 
below about O.SMq. From extrapolation of simple an- 
alytic models for the surface density distribution, we 
infer that about half of the lower mass cluster members 
lie outside our surveyed area, whilst the vast majority 
of higher mass stars are included. If this is confirmed 
by wider surveys for cluster members, then the whole- 
cluster mass and luminosity functions of NGC 2516 
and the Pleiades could be very similar at least down 
to 0.3Mq. The cluster core radius increases towards 
lower masses at a faster rate than if the segregation 
were simply due to equipartition. This may indicate 
some remnant primordial mass segregation in NGC 
2516, where the high-mass stars were initially more 
centrally concentrated than the low-mass stars. 

4. The binary fraction of A to M-type stars, with mass 
ratios of 0.6-1, is 26 ± 5 percent in NGC 2516, which 
is comparable with the same statistic for the Pleiades 
and field stars. This is a lower limit to the total bi- 
nary fraction, which could be as high as 65 to 85 per- 
cent, depending on what form the mass ratio distribu- 
tion takes. Unresolved binarity has only been partially 
taken onto account in our mass function determina- 
tions, in the sense that we have identified those systems 
composed of two stars of roughly equal mass and have 
not treated them as a single star with slightly higher 
mass. Mass functions determined in this way at least 
have the merit of comparability with estimates derived 
from photometric data on other clusters in the litera- 
ture, but we recognize that the slopes of the true stellar 
mass functions may be a little steeper (a increases by 
^ 0.05 — 0.1 for high mass stars and r^ 0.3 — 0.4 for 
lower mass stars), because low mass stars could still 
be hidden in many systems. 

5. The total mass of the cluster for stars with M > 
0.3Mq is at least 950Mq and probably as high as 
1200Mq, depending on the metallicity, binary fraction 
and mass-ratio distribution. Correcting for the frac- 
tion of stars that lie outside our survey, but within the 
likely cluster tidal radius, may increase this estimate 
to around 124O-156OM0. NGC 2516 is therefore about 
twice as massive as the Pleiades. 

6. If the whole-cluster mass functions of NGC 2516 and 
the Pleiades are similar, then we expect about 360-440 
brown dwarfs in NGC 2516. Extrapolation of the mass 
functions derived from our data suggest that 100-220 
would be in the area surveyed in this paper. 
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